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CN , Using a quantum theory for an ensemble of three-level atoms (A) placed in 

D . an optical cavity and driven by electromagnetic fields, we show that the long- 

. lived spin associated with the ground state sublevels can be squeezed. Two 

. kinds of squeezing are obtained: self-spin squeezing, when the input fields 

. are coherent states and the atomic ensemble exhibits a large non-linearity; 

o ; 

On , squeezing transfer, when one of the incoming fields is squeezed. 

(N ■ 

o : 

■ 

C ■ I- INTRODUCTION 

a • 

^ : 

^ ' In this paper, we present calculations performed on a model system consisting of three- 

'X' 

' level atoms that interact with two fields. The atoms are placed in an optical cavity which 
ensures a sizeable level of interaction, in contrast with single pass schemes. Under such 
conditions, the validity of quantum fluctuation calculations based on quantum Langevin 
equations is well established. We will consider a set of three-level atoms in the so-called A 
configuration, consisting of two ground state sublevels and one excited state. 

It is well-known that, in the case of large one-photon detuning (Raman scheme), terms 
involving the ground - excited state coherence and excited state population can be adiabat- 
ically eliminated. The three-level system is then equivalent to an effective two-level system 
in which the spin associated with the ground state sublevels may be squeezed. Because of 
its simplicity the effective two-level system provides a good understanding of the physical 



phenomena responsible for the creation of atomic squeezing and allows one to carry out 
analytical calculations of optimal squeezing. 

We study two regimes. In the first, the atoms interact with coherent light. One can 
achieve atomic self-spin squeezing if the non-linearity of the medium is sufficiently high. This 
atomic squeezing originates from the same physical process that gives rise to the squeezing of 
the field exiting the cavity, predicted theoretically in Refs. ||l|,@J^ and observed experimen- 
tally a little later In the early 1980's, it was conjectured that atomic spin squeezing 
appears as a counterpart of squeezing of the electromagnetic field . However, the quan- 
tum noise reduction on atomic variables computed in several papers can be obtained by a 
rotation of the atomic variables of a two-level system interacting with a coherent field. As 
shown in in order to be relevant, spin squeezing should be computed in the plane 
orthogonal to the direction of the mean spin. We have performed full quantum calculations 
in the relevant basis. Here we generalize the results obtained previously on two-level systems 



and three-level (V) |T^ systems to the A system. 

In the second regime, the atoms interact with a squeezed incoming field. This scheme 
was proposed recently and experimentally tested in a single pass configuration to produce 
spin squeezing between two excited levels [0,^. In |T^, it was shown that the cavity 
configuration should allow one to enhance the squeezing transfer between two-level atoms 
and field. In both cases the transfer was limited by the vacuum field fiuctuations. We find 
here that the A scheme allows for a quasi-perfect squeezing transfer from the field to the 
atoms. Unlike configurations in which the squeezing occurs between two excited states or 
between a ground and excited state, the squeezing is created between two long-lived ground 
states. Such a system offers significant advantages from an experimental perspective. 
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II. EQUATIONS FOR ATOMIC FLUCTUATIONS 



A. Three-level atoms interacting with two fields 

We consider an ensemble of three-level atoms, the three levels being labelled 1,2,3 in a A 
configuration (Fig. 1). The atoms interact with two light fields: an intense classical pump 
field Ai, having frequency cui and a quantum probe field A2, having frequency 002- Field Ai 
transverses the medium in a single pass and drives the 1 — > 3 transition, while field A2 is 
confined to an optical cavity and drives the 2 — 3 transition. The detunings from atomic 
resonance are Ai = cusi — cui for the pump, and A2 = CU32 — ^2 for the probe. The cavity 
resonance frequency closest to the probe frequency is lOc- The cavity detuning for the probe 
field is Ac = ci;2 — tUc- The incoming quantum field that is coupled into the cavity is denoted 
by A^^2 ■ The intensity of field Ax is assumed to be much greater than that of the quantum 
field. 

The three-level system is described using nine collective operators for the atoms of 
the ensemble, the populations of levels |1) , |2) et |3), 

N N N 

Hi = E , = E |2), (21, , Ha = E |3). (3|, , (1) 

i=l i=l i=l 

the components of the optical dipoles in frames rotating at the laser frequencies, 

N N 

^iW = E|l).(3Le^'^^ AW=E|2).(3Le^'^^* (2) 

i=l i=l 

(and their hermitian conjugates), and the operator associated with the coherence between 
levels |1) and |2), 

N 

m-J:\^)^{Mie'^'''-^^^' (3) 

i=l 

(and its hermitian conjugate). 

The Hamiltonian for the atom-field system is given by 

H = h(jj2i'si2 + h(jj3iiis + huj2a\a2 + hg{P^A2 + AIP2) + h{niPi + niPl), 
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where g — Sod23/h is a coupling constant, assumed real, associated with the (quantum) 



cavity field, Sq — y ^c<j2/2eo»Sc is an amplitude appearing in the equation for the cavity 
field operator E2 = Sq (e'^'^^t^^ + e'^^^Al^ A2 and A2 are field creation and annihilation 
operators of having commutator 



K is the cavity decay rate, t = L/c, where L is the cavity length, S is the cross sectional 
area of the cavity field, Qi = Eidis/h is a couphng constant associated with the classical 
pump field, Ei{t) = (e~"^i*£Ji + e'^^'^^E*) , ^23 and dis are dipole moment matrix elements. 
Note that f^i has units of frequency, while g and A2 have units of (frequency) The decay 
constants associated with dipole operators Pi and P2 are taken to be equal and denoted by 
7. In order to account for the finite lifetime of the two ground state sublevels 1 and 2, we 
include in the model another decay rate 70, which is much smaller than 7. For example, 
7o can refiect transit-time broadening in the system; as such, it affects level 3 as well as 
levels 1 and 2. Typically, the atoms fall out of the interaction area with the light beam 
in a time of the order of tens of milliseconds, whereas 7 is of the order of the MHz for 
the excited state. We also include incoherent pumping of levels 1 and 2 with rates Ai and 
A2, respectively, to allow for ground state atoms entering the interaction volume. We will 
neglect any fiuctuations in the total number of atoms in the interaction volume and assume 
that Hr = Hi + n2 + Ha = (Ai + A2) /70 = N. 

The time evolution of atomic polarizations and populations is obtained by adding to the 
Heisenberg equations of motion, ihO = [0,H], terms corresponding to the Langevin forces 
associated with the decay of these quantities, 
dUi 



= iO*Pi - i^iPl + -fUs - 7oni + Ai + Fn, (4) 
igAlP2 - igA2Pl + 7n3 - 7on2 + A2 + Fn^ (5) 



dt 

dTl2 
dt 

^ = -{1Q.IP1 - tQiPl) - {igAP2 - igA2Ph - 27n3 - 7on3 + (6) 



dP 

' -(^Ai + 7)Pi + tn^iU, - Us) + igA^P} + Fp, (7) 



dt 

dP2 ^ 
dt 

dPr 



dt 



{iA2 + i)P2 + igA2{Ti2 - Hs) + l^iPr + Fp, (8) 
= - (70 - i5) Pr + tniP2 - igA2Pl + Fp^ (9) 



where 



5 = Ai - A2 

is the detuning between the ground state sublevels. " In" terms for resuhing from sponta- 
neous emission have not been included in Eq. (|^), owing either to polarization selection rules 
or a hyperfine separation of levels 1 and 2 that is much greater than 7. To these equations 
must be added the evolution equation for the cavity field. 

Changes of the intracavity field result from the incoming field A™, losses through the 
coupling mirror, and the source field associated with the atomic polarization. The evolution 
equation for the intracavity field is W2\ 



^ = -(/. + .A.) A2 + ^P, + .P^r (10) 
at T \ T 

As is evident from this equation, fluctuations of the incoming field give rise to a Langevin 
force for the cavity field. 

Our aim is to obtain the fluctuations of the spin operators associated with levels 1 and 
2. For this, we need to determine the fluctuations of operators P^, P}., IIi and 112. To obtain 
equations for the fluctuations, we linearize equations (§-1) , and (|^^, |10|) and their hermitian 
conjugates. We define a 10- dimensional vector |^(t)] as 

m] = \Mt), 4{t), Pl{t), Pht), P2{t), PKt), Pr{t), P}{t), S^^{t), S,2{t)f 

where 

5,1 = (Hi - Hg) /2; S,2 = (n2 - Hs) /2 
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and then set 

\m] = mm + mt)] 

The steady state values, {^(t)) , are obtained from Eqs. (0) in steady state, using 

the fact that the mean value of the Langevin operators are zero; however, these analytical 
solutions are of marginal use. 

The fluctuation vector obeys an equation of motion 

dmt)] 



dt 



[B] mt)] + \m] (11) 



where [B] is the linearized evolution matrix of the atom-field system and the column vector 
contains the Langevin forces: 



\F^{t)] = \pKlr6A^\t), ^2K/T64^\t),Fp,{t),Fp,{t), (12) 
Fp, (t) , Fpt (t) , Fp^ (t) , Fp, (t) , Fs^, (t) , Fs^, {t)f 

The correlation matrix [G(t)] of the fluctuations is: 

m)] = mt)]mo)\ (13) 

The variances are obtained from the zero time correlation functions, contained in the matrix 
[G(0)] which satisfies: 

[B] [cm + [cm [B]' = [D] (14) 
where [D] is the correlation matrix of the Langevin forces : 

{\m]mt')\) = [D] s{t-t'). (15) 

The sub-matrix [GclO)] is defined as the 4x4 lower diagonal block of [G'(O)] and contains 
the variances of Pr, Pj., S^i, and 5*22 ■ 
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B. Spin squeezing 



1. Definition 

In the same way as a squeezed state of the electromagnetic field is defined by comparison 
with a coherent state, a squeezed spin state will be defined as having fluctuations in one 



component lower than that of a coherent spin state [|TT| . A coherent spin state for atoms 
is deflned as a product state of N uncorrelated spins, in which the state vector for the i-th 
spin is an eigenstate the individual spin operator in the {6, 0) direction: 

o'e,(t> i = ^xi sin 9 cos (p + ayi sin 6' sin + a^i cos 6 (16) 

having eigenvalue +1/2, with i = (cr+ i + crl j)/2, Cy i = j — crl i)/2i. The coherent 
spin state can be obtained as a rotation from the state having all spins aligned along the 
z axis. This coherent spin state is an eigenvalue of the collective spin operator Se^,/, = 
Y.i=i,N ^e,<j> i, with eigenvalue S = N/2 ||18[. It satisfles the minimum uncertainty relationship 
with fluctuations equally distributed over any two orthogonal components normal to the 
{0,(f)) direction, the variance of which is equal to S/2 = N/A . If one can squeeze the 
fluctuations of the total spin within the plane orthogonal to the mean value, it will result in 



noise reduction in spin measurements. The condition for spin squeezing is then |T9 



< 1(5^)1/2 (17) 

where the axes have been rotated in such a way that the Z axis is in the direction of the 
mean spin and a represents a direction in the X, Y plane. The quantity (Sz) is then the 
mean value of the spin and Sx and Sy have zero mean values. Spin squeezing can occur 
only for a spin ensemble with > 1 since it implies the emergence of quantum correlations 
within the spin ensemble. 
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2. Minimum variance calculation 



We calculate the variances ASx and ASy of the spin variables in the new reference 
frame. For this, we perform a rotation defined by angles (p around the z axis and 9 around 
the Y axis (defined by the first rotation) such that 



cos 61 = (S^)/ I {S) I COS0 = {S^)/S^ 



(18) 



with I {S) 1= ^{Sx^ + {Syy + (5*2)2 and = ^ {SxY + (5'j/)^, where the mean values are 
the solutions of the steady state equations. A spin component in the X, Y plane making an 
angle a with the X axis has a variance given by: 



AS a = cos^ a ASx + sin^ a ASy + sin 2a Re {ASxy) 



(19) 



with AS"^,^ =< SSn{0) SSu{0) > {iJ,,u — X, Y). As a consequence the values ckq of a for the 
spin components having maximal and minimal variances satisfy: 

2 Re (ASxy) 



tan 2ao 



(20) 



ASx - ASy 

In order to investigate squeezing, we compare the minimal variance to K^'^)] /2. The cor- 
responding normalized variance thus obtained is called AjSmin- Spin squeezing is achieved 
when A5'mm < 1- 

For the three-level system, the spin fluctuations corresponding to the atomic coherence 
between levels 1 and 2 are given by lower diagonal block [(^0(0)] of the correlation matrix. 
To go into the relevant basis, one must first perform the transformation from Pr, Pj, S^i , 

Sz2 to Sex, Scy, Scz with: 

n2 - Hi 



P _L Pt p — pi 

^ r \ r q _ ^ r ^ r n 



For this we define matrix : 



2i 



'^1100^ 



-i i 
-2 2y 



(21) 



(22) 
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Then one goes into the basis Scx, Scy, where the mean spin direction is along OZ, by using 
the [-R2] matrix: 



V 



[R2 



with angles 9 ei ip given by Eq.(|T 
The atomic correlation matrix 



^ cos 9 cos cos 9 sin ip — sin 9 ^ 



sm ip cos ip 







(0) in the Scx, Scy basis is 



Gj{0) = [R2] m [G,(0)] [R 



The variances in the XY plane are then 



AS, 



cX 



1,1 



AS, 



cY 



cm 



2,2 



cXY 



cm 



1,2 



The minimal variance ASc min corresponds to the angle given by Eq.f^Ol). 



III. ADIABATIC ELIMINATION IN THE RAMAN SCHEME 



A. Effective two-level system 



(23) 



(24) 



(25) 



In the limit that | Ai^2| ^ 7? one can adiabatically eliminate the ground-excited state 
polarizations. In this limit, the excited state population is much less than unity and can be 
neglected. We assume also that \6\ = \Ai — A2I ^ |A| = (Ai + A2)/2, since the values of 
1^1 to be considered are of order of 70 ^ |A|. Performing this adiabatical elimination and 
neglecting terms of order 7/A^, one obtains simplified equations for the sublevels variables 
S+ = Pr, S. = P} , S, = (n2 - ni)/2 and n = n2 + Hi = iV: 

dS4 



dt 



dS- 
~dt 



-{^o~^^)S+ + 2i~gA2S, + Fs^ 



-{lo + i5)S^-2i~gA\S, + Fs 



(26) 



(27) 



(28) 



^ = -7on + A (29) 

where A = Ai + A2, 

^ = ^ + ^ -K- 

is the effective atomic detuning corrected for hght-shifts, and 

9 = 9^ 

is an effective couphng constant. The evolution equation ( p!OD for the quantum field becomes 

^ = -(k + tAM2 + i^S^ + (30) 

In principle, the atomic fluctuation operators contain terms related to repopulation of the 
ground states via spontaneous emission; however these contributions are down by (7/A)^ 
and are neglected. As a consequence, the fluctuation operators are associated only with the 
decay rate 70 and the repopulating processes; we shall see that the diffusion matrix formed 
from these operators has a fairly simple form. 

The problem has been reduced to that of an ensemble of two-level atoms interacting with 
one field in an optical cavity whose effective Rabi frequency is fie// = 9 (^2) = 9^ (^2) ■ 
We will use these simplified equations for the spin squeezing calculations and give a physical 
interpretation to our results. Note that, beside the incoherent pumping terms, our two-level 
system differs slightly from a classical "excited-ground state" two-level system. The decay 
terms in (|2^-|2^^Bp are all 70; the coherence relaxes at the same rate as the populations. 



This remark, as simple as it may seem, critically affects the steady states that are allowed 
by the system (see Sec. IV C). Note also that the atomic noise scale is no longer 7, but 
7o, hence much smaller. We thus expect the atomic noise not to destroy the squeezing too 
much. 
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B. Steady state and bistability 



The steady state values are obtained solving Eq. P6| - p7| - p8l - pO| ) in steady state, using the 
fact that the mean value of the Langevin operators are zero. Explicitly, one finds 



(S-) 



A2 -Ai 
7o 



P2 



6 



1 + 52 + 4J2 



A2-A1 
7o 



/?2 



6 



1 + 52 + 4/2 



A2 - Ai 1 + 6' 



270 1 + 5' + 4/2 
with the dimensionless variables (in units of 70) 



(31) 



(32) 



(33) 



7o 



132 



7o 



Pi 



(34) 



For certain values of the parameters, the system exhibits bistability. From ( pOD in steady 
state, one derives the mean value of the intracavity intensity 

2C V ( . 2C5 \ 



Tin 
-'2 



+ Uc + 



(35) 



1 + 52 + 4/2; V l + S' + AhJ 

where Sc = Ac/ n is the dimensionless cavity detuning, and C is an effective cooperativity or 
strong-coupling parameter 

g' A2 - Ai 



C 



(36) 



2Kr7o 7o 

We will limit our discussion to an incoherent pumping scheme in which Ai = and A2 = iV7o, 
corresponding to absorption for the probe field. We will show in the next sections that it is 
indeed the most favorable case for squeezing. The effective cooperativity is then related to 
the usual cooperativity parameter C = g'N/2K,T'y by 



C = C 



A 



1_ 

7o 



(37) 



The cooperativity parameter determines the efficiency of the non-linear effect. If C - or C 
- and are large enough, the intracavity field exhibits a bistable behavior. It has been 
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shown that the field quantum fluctuations are maximal around the lower turning point of 
the bistability curve, allowing for the outgoing field to be squeezed. In Ref. [0, it has been 
demonstrated that this point is also interesting for atomic squeezing. Therefore, we will 
study the fiuctuations in the vicinity of this turning point. 



C. Linearization and difTusion matrix 



The linearized equations obtained from (p6| - pOD are 



"^^^^ -(70 - ^~S)5S^ + 2t~g{S,)6A2 + 2i~g{A2)5S, + F, 



dt 



dSS_ 
dt 



-(70 + ^S)SS. - 2t~g{S,)6Al - li-g^A^YbS, + Fs 



dSSz 
dt 



-loSS, - i~g{{S-)SA2 - {S+)6Al - (A2)*(55+ + {A2)6S^) + Fs 



(38) 
(39) 
(40) 



d6A2 
dt 



[k + iA^)6A2 + ihs+ + \—6A^:^ 

T \ T 



We can write Eqs. (|38| - ^TD and the hermitian conjugate of (|4l|) in matrix form as 

d mt)] _ 



where 



dt 



[B] mt)] + \F^] 



[B] 



mt)] = [6A2it),6AUt),5S4t),5S4t),5S,it) 

^ {K + iA^) -ig/r 

(k-zAc) ig/r 

-2tg{S,) ilo-^5) -2tg{A2) 

2ig{S,) (70 + ^5) 2t~g{A2)* 

\ igiS-) -ig{S+) -i9{A2)* ig{A2) 70 



(41) 



(42) 



(43) 



and 
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FsAt), Fs_ (i), Fs. (i) 



T 



(44) 



The covariance matrix is defined by 



[G{t)] = mt)] mo)\ 



(45) 



and the diffusion matrix by 



{\m]m')\) = [D] 6it-t') 



(46) 



The values of the atomic diffusion coefficients can be derived from the quantum regression 
theorem or operator identities. The complete diffusion matrix is given in the Appendix. 
The variances of the spin components and their correlation functions are the elements of the 
zero time correlation matrix [G(0)], which satisfies [§ 



Inverting Eq. (0), one obtains [G(0)] , and, consequently, the spin variances. We proceed 
with the calculation of the minimal variance as in the first Section. 



We first study the case of a coherent input field A^^. We have studied a wide range of 
parameters. First, in Fig. 2, we have plotted the minimum variance AS'^in as a function of 
the cooperativity C for a typical steady state point (5 = 10, 6c = 0, I2 = 25.2). 

We see that the squeezing naturally increases with the cooperativity, the non-linear 
interaction increasing, but saturates a little under 30 % when C exceeds 100. As it is 
a relatively accessible value for experiments, we will often choose this value for C in the 
following. 



[B] [cm + [cm [B]' = [D] 



(47) 



IV. SELF SPIN SQUEEZING 



A. Minimal variance calculation 
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Then, for a fixed cooperativity (7, we fiave optimized, for each value of 5, the values of Sc 
and I2 yielding the best atomic squeezing, that is the lowest A^min- The results are reported 
in Table I for (5 = 100. 



s 





5 


10 


15 


20 


Sc 





-0.2 





-0.4 


-0.2 


h 


0.25 


6.5 


25.2 


56.5 


100 




0.713 


0.716 


0.72 


0.72 


0.728 



Table I. Minimal spin variances for the effective two-level system as a function of 5 (C 
is equal to 100, the other parameters are adjusted for optimal noise reduction). 

We see that the squeezing does not vary much with 5, and is always close to 28-29%. 
Increasing the cooperativity does not improve the squeezing much: one saturates at a squeez- 
ing value of almost 30%. The fact that the minimum variance is rather independent from the 
triplet ((5, (5c, I2) when the cooperativity increases shows that there exists an optimal turning 
point for the system that reaches a steady state independent from the above mentioned 
triplet as soon as the cooperativity is large enough. Furthermore, the minimal variance 
seems to be rather indifferent to the atomic noise for large values of C. The following 
sections aim to explain these a priori surprising results. 



B. Outgoing field spectrum 

As mentioned above, with the existence of atomic squeezing, we can expect squeezing 
for the outgoing field A^^. To analyze field squeezing, a good variable is the noise spectrum 
in the frequency domain obtained from the variance matrix \yout {'^)] defined by 



(48) 



where 



(49) 
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is the vector Fourier transform of I^Ag"* (t)] . The minimal and maximal component of the 
noise spectrum are then 

Sout min (t^) = [Vout + [Vout (t^)]2,2 ^ 2 \[Vout (t^)]li,2 (50) 

Sout max 

(cu) = [Vout 1 + [Vout {^)]2 2 + 2 \ [Vout 2 (51) 

We have plotted in Fig. 3 these two components for a given set of parameters of Table I. 
As expected, the field is squeezed (maximum squeezing around 12 %), but the squeezing 
occurs over a frequency range having width IOO70, much larger than the value of order 70 



one might have expected from Eq. (|38| - pD where the atomic time constant is l/7o- This is 
a consequence of the field-atom interaction in the "bad cavity limit", as will be explained 
in the next section. 



C. Adiabatic elimination of the intracavity field 

Unlike the three-level system, the effective two-level system we consider is in what is 
called the "bad cavity limit": k ^ 70, the field evolves much faster than the atoms. We 



can therefore adiabatically eliminate the field operators in the linearized Eqs. ( p8[]4l|) . One 
obtains a new set of equations 

d\6S] 



dt 



with \SS] = [SS+,6S_,SS,\^, \F^] 



-70 [53] \6S] + \F's] 



pi pi pi 



[Bs 



1-16 + AC 







l+i6 + 4Cj^^ 2t(32 



2C 



l+i<5c 



where 



(52) 



(53) 
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In these equations there appears a new decay constant 

Sz 



7 



7o 



1 + 4C: 



1 + 61 

which is approximately equal to C70 if C ^ 1. If we take the values of the parameters of Fig. 
3, we find 7' = IOO70 which corresponds exactly to the frequency bandwidth over which the 
outgoing field is squeezed. The atomic spectrum is broadened as a result of the interaction 
with the field. The new diffusion matrix defined by {\F'{t)] [F'{t')\) = [D'] 6{t — t') is given 
by 



[D'] = 7oiV 



[Dat] + C 



V 



V 



(54) 



Asl -2s,s^ 

—2szS- 

where [Dat] is the atomic diffusion matrix calculated from the quantum regression theorem 
(see Appendix). 

Equation (|54D shows that [D'] is the sum of the usual atomic diffusion matrix '-j^NDat, 
which is proportional to the atomic decay rate 70 and the number of atoms A^, and another 
matrix, which represents the contribution of the field fluctuations, proportional to A^^, since 
C (X N . For large values of the cooperativity, the main contribution to the noise originates 
from the field. This conclusion is valid only if 



7o < C70 < k; 



(55) 



if C becomes too large, the atoms evolve as fast as the field, and the adiabatic elimination is 
no longer justified. When inequalities (^) hold, one is able to perform analytical calculations 
for the minimum variance, considering only the larger terms in ( ^4]) and (|53|), and using an 
optimal bistable point as in Table I. It is possible to show that the lowest minimum variance 
^^min that can be obtained is = 0.707, in good agreement with the results announced 

in Sec. |1V A| . The optimal bistable point is obtained considering the most favorable steady 
state for squeezing. Since the squeezing originates from atom correlations we expect it to be 
maximum when the modulus of coherence 5*+ is a maximum. From the steady state equations 
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( pl| - |3^|33D , it is straightforward to show that the maximum value of |S'+| is when 
4/2 = 1 + 5^. In this case, Sz = N/A and, integrating (p2D, one finds, for C 1, 

A^^in = ^ + 0(l) (56) 



We see that the self squeezing limit in this system is approximately 30 %, as was found in 
the numerical simulations. 

D. Noise contributions to the atomic variance 

To confirm our statement that the squeezing can be traced to fluctuations in the cavity 
field, we have computed the different contributions to the atomic noise. Indeed, instead 
of calculating directly the atomic variances from the time domain equations, it is possible 
to derive the atomic noise spectra in the Fourier domain and then integrate the various 
contributions over the frequency to get the variances. The linear response theory allows for 
calculating the response of the atoms to the field excitation. The fluctuations of the incident 
field cause the intracavity field to fluctuate, which induces fluctuations for the atomic vari- 
ables. These fluctuations are coupled back to the field by polarization fluctuations, hence 
the non-linear coupling between the intracavity field 15^42] and the spin fluctuations 155*]. 
To these sources of fluctuations, we must add the coupling with the vacuum surrounding 
modes of the cavity and the noise due to the decay and pumping of the sublevels. Taking 
all these fluctuations into account, the atomic variance matrix [Vaj [u)] is the sum of four 



contributions [12 



[Vat (^)] = [Vf {UO)] + [K {UO)] + [Vd^p {uo)] + [Vnt (^)] (57) 

in which 

1. \Vf {uj)] represents the contribution of the incident field whose fluctuations are modified 
by the atoms 

2. (a;)] is the contribution of the vacuum fluctuations 
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3. [Vdip (cu)] corresponds to the fluorescence emitted in the cavity mode 



4. [Vint (i^)] is the interference term between the vacuum fluctuations and tlie dipole 
fluctuations. 

In Fig. 4, are plotted the various contributions to the minimal component of the spin 
noise spectrum. It is clear that the atomic noise is due entirely to the fluctuations of the 
incident field, whereas the sum of the contributions from the other terms is negligible. We 
have integrated the different spectra to yield the contributions to the variances ASmin and 
present the results in Table II. 





A^, 


AS dip 


ASint 


A ■ 


ASf/AS^i^ 


0.701 


1.414 


1.373 


-2.772 


0.716 


97.9 % 



Table II. Contributions to the atomic variance (value of the parameters: C = 100, 6 = 12, 
4 = -0.2, I2 = 40, 7o = 7/1000, K = 27 = 5.2MHz) 

We see that the contribution of the field represents about 98 % of the total variance for 
the minimal component, whereas all the other contributions amount to only 2% of the global 
noise. Consequently, the A scheme provides a non-linear regime in which the field imprints its 
fluctuations on the atoms and in which all the other causes of noise are negligible compared 
to the non-linear interaction. We can thus easily reach the squeezing limit predicted in the 



previous section, for much smaller values of the cooperativity (50-100). In [13|, for a classical 
ground-excited state two-level system, the self-squeezing limit was 50 %, but, in order to 
approach that limit, much larger values of the cooperativity were required (1000 to 10000). 



V. SQUEEZING TRANSFER 

The second configuration that we have studied is the case where the input light is a 
broadband squeezed vacuum. In this case, (^2") — 0- Note that we are no longer in a 
non- linear regime as in the previous section. We have shown that, for the 3-level Raman 



scheme considered in this paper and for the optimized transfer conditions given in Ref. ||16|| , 
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the squeezing transfer from the field to the atoms is almost 100 %, even though we are not 



in the strong coupling regime as in [0. The only approximation made is the same as in 
Ref. |T6|: the intracavity intensity stays negligible and the steady state mean spin is then 



{Sz) = N/2 ("small angle approximation" |TT1). Of course, even though {A^^) = 0, the 
average number of intracavity photons is not 0, but sinh^r for a squeezed vacuum, where r 
is a squeezing parameter. If the average photon number is much less than A^, however, one 
can set 

(^+) = (^-) = 0, {S.) = ^. (58) 

Under these assumptions, the calculation of the minimum variance can be done analytically 
T^, and one can deduce the optimal transfer condition 



5 = 4 = 0. (59) 

In this case, the minimal variance can be expressed as a function of the effective cooperativity 
C, the ratio p = 7o/k and the squeezing rate e"^'' as 

9C 

A5'^in = l ^(l-e-'n. (60) 

il + p)il + 2Cy ^ ^ ' 

As we pointed out, we are not in a strong coupling regime, which corresponds to pC ^ 1. In 

our case, C ~ 100 and p = 7o//« ~ 1/2000 ^ 7/^ (for a ground-excited state (g-e) system 

with K = 27, p = 1/2 and pC ^ 1). For large values of C, the factor in front of (1 — e"^'') 

is close to unity. If we define the transfer efficiency 77 by 

1 — A '? ■ 

we find, for large C, 

1 - (62) 



(l + p)(l + ^) (1 + P)V 2CJ- 

In particular, for infinite cooperativity, we find the same limit of maximal squeezing as in 
|[T6| , but in a regime for which 70 <^ k (p <^ 1, ~ 1 — p). In Fig. 5, we have plotted 
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the minimum variance versus the squeezing of the input field Rin = 1 — e"^*" , for C = 100 
and two different values of p corresponding to the case of our effective two-level system 
(p = 1/2000), and the case of an g-e system (p = 1/2). In particular, the A system enables 
one to reach much higher squeezing values than in a classical g-e two-level scheme. To 
understand the influence of the cooperativity on the transfer, the efficiency versus C is 
plotted in Fig. 6 in the two cases discussed previously. The flrst conclusion is that the 
transfer is close to 100 % in the regime considered (p = 1/2000) and is not limited by the 
vacuum noise as in a g-e transition [T^. Second, the efficiency quickly saturates with C; a 
cooperativity of 100 enables a transfer of 99.5 % of the squeezing from the field to the atoms 
(p = 1/2000). We find again that the transfer efficiency is excellent for very reasonable 
values of the cooperativity (~ 100), which is usually the limiting experimental factor. 



VI. VALIDITY OF THE TWO-LEVEL MODEL 

In this section, we discuss the validity conditions for our effective two-level model and 
see how the model can be modified slightly to extend its range of validity. In carrying out 
the adiabatic elimination that took us from the full, three-level equations to the two-level 
equations, two basic assumptions were made. First, it was assumed that optical pumping 
terms of order 

could be neglected. Second, it was assumed implicitly that the probe field was sufficiently 
weak to neglect terms of order 1^2/^1 1 , where 

If one adiabatically eliminates both the excited state - ground state coherence and the 
excited state population operators from Eqs. (|-|^), the resulting contributions to the time 
evolution of the ground state operators of optical pumping, transit time decay, and incoher- 
ent pumping is 
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S^ = - (70 + Tp) + NTp/2 + (A2 - Ai)/2 (63a) 
S+ = - (70 + Tp) 5+ - ^p^]2/^^l (63b) 

where terms of order Fp 1^2/^1!^ have been neglected. Aside from the term of order 
Fp 1^2/^1! J which we neglect for the moment but return to shortly, these equations have 
basically the same structure that we used for our effective two-state model. If Fp <^ 70, the 
two models coincide. For arbitrary ratios of Fp/70, one can still use an effective two-level 
model if the assignments 

7o (two-level) = 70 (three-level) -|- Fp (64a) 
Ai (two-level) = Ai (three- level) (64b) 
A2 (two-level) = A2 (three- level) + A^Fp (64c) 

are made [recall that Ai -|- A2 = A^7o]- To check this hypothesis, we set 70 (three- level) = 
Fp = 7/1000, so that Fp is no longer negligible compared to 70, and we compare the two- and 
three-level calculations for self-squeezing when A = IOO7, fii = v^T07. For these parameters, 
C = C = 100 and 1^2/^1! is typically of order 0.01 at the point of optimal squeezing. The 
results, displayed in Fig. 7, are in good agreement for both models. We have generalized 



the results obtained in Sec. |IV A| for Fp/70 <^ 1 to arbitrary ratios Fp/70. We would like to 



point out that the steady state values (0), (|3^) , ( ^3|) are modified by the presence of Fp, 



and the relation 4/2 = (1 + Fp/70)^ + 5^ must be satisfied in order to have the same optimal 
value of l^+l = N/A. 

To achieve better agreement between the effective two-level and full three-level calcula- 
tions, it may also necessary to include the |f22/fii| term in Eq. (|63b|) . In the previous cases, 
this term can be neglected, but, for closed systems [Fp ^ 70 (three-level)], the term linear 
in 1^72/^1! in Eq- (63b) can slightly increase the value of KS*-!-)! above the maximum value 



of A^/4 it has for open systems [Fp ^ 70 (three- level)]. As such the value of ASmin can be 
reduced below 1/V2- This feature is seen in Fig. 8, where inclusion of the linear term in 
a corrected two-level model brings the results into good agreement with the full three-level 

21 



calculation. Squeezing values of 35% can be reached in such a regime. The corrected two- 
level model also includes a slight correction arising from a term linear in 1^2/^1! that we 
neglected in going from Eq. (|T0|) to ( pO| ) for the field evolution. When this linear term is 
included, Eq. ( |30| ) is modified as 



■ A \ A '''9 



dt ' ' ' ' T 



K + 2 A,) A2 + - 5+ + n^/n, n2 +\ — A*2" (65) 



2k 



r 



The simplified two-level model considered in the earlier sections provides a good under- 
standing of squeezing in a three-level system and allows one to optimize the parameters in 
a very simple way. It can be brought into full agreement with the three-level model when 
1^^1,21 ^ 7) l^il and |i72/^i| ^ 1) provided one uses the prescription (|6^) for the decay and 
incoherent pumping and Eq. (|65|) for the field evolution. Finally, we would like to point out 
that, although the value of 70 itself is not critical, the ratio of the decay constants of the 
coherence S± to the population difference is important, since it determines the maximal 
coherence that can be obtained from the system, hence the maximal squeezing. Changing 
this ratio from 1 (as in our simple two-level system) to 1/2 (as in a classical e-g system) 
would enable one to achieve a 50% squeezing limit. 



VII. CONCLUSION 

Using a quantum model for an ensemble of three-level atoms in a Raman (A) type 
configuration, we have derived the atomic spin fluctuation spectra and variances and have 
shown rigorously the occurrence of spin squeezing between the two ground state sublevels. 
Unlike configurations in which the squeezing occurs between two excited states or between 
a ground state and an excited state, the squeezing obtained here has a long lifetime, does 
not require a very strong coupling and, as a consequence, presents substantial advantages 
for experimental realization. 

Spin squeezing may occur in two different regimes. In the first one, the non-linearity 
of the atomic ensemble is exploited to squeeze the intracavity field, which in turn imprints 
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squeezing on the atomic ensemble, yielding self-spin squeezing. In the second one, the atomic 
ensemble has a linear behavior. It cannot create squeezing in the intracavity field. However, 
if one of the incoming fields is squeezed, the atom-field coupling in the cavity enables a 
transfer of squeezing of almost 100 %. 
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IX. APPENDIX 



We give here the expression of the diffusion matrix appearing in (|^) for the two-level 
system. The matrix elements of the upper 2x2 diagonal sub-block of [D] are the correlation 
functions of a broadband squeezed field, equal to the one of a single mode field squeezed by 



a factor e ||T6| 



< >= cosh2(r) S{t - t') 



(66) 



< M*"(t)M™(t) >= l/2sinh(r) e''5{t - t') 



(67) 



< 6A'''\t)6A'"\t) >= l/2sinh(r) e-'^6{t-t') 



< SA'''\t)6A'''{t) >= sinh'(r) 5{t - t') 



(68) 
(69) 



leading to a field diffusion matrix 



r 1 2/? 



cosh^(r) l/2sinh(r) e *^ 
1/2 sinh(r) e*^ sinh^(r) 



(70) 
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The matrix elements of the lower 3x3 diagonal sub-block of [D] are the correlation functions 
of the atomic noise operators appearing in Eqs. (p^)-PO|). They were evaluated with the 



Einstein generalized relations ||T^. and grouped in the atomic diffusion matrix [D, 



at 



at 



iV7o 



1 _|_ Ai— A2 







:i + Ai-A2 







;i + Ai-A2)4^ 



(71) 



\2; — 

+ (-1 + Ai-A2)^ 
2 v-l + Ai-Aa)^ | + (Ai-A2)s, 
where Ai,2 = are the dimensionless incoherent pumping terms for one atom and 

s±,z = {S±,z) /N are the steady state values for one atom. The other elements of [D] are 
equal to zero since there is no correlation between atomic and field fluctuations at the same 
time. Therefore, 



[D] 



( [Dj] ^ 



V 



[Dat] 

where [Dat] and [Df\ are defined above (|70])-(^). 



(72) 



Fig 1. Three-level system in a A configuration. 

Fig 2. Minimal variance versus effective cooperativity C for a typical steady state point 
[5 = 10, 5c = 0, /2 = 25.2). 

Fig 3. Minimal (below) and maximal (above) spectra for the outgoing field for the same 
bistable point as in Fig. 1 (the frequency unit is 70). 

Fig 4. Contributions to the atomic noise spectrum (black): the field contribution (light 
grey) is predominant over the sum of the three other contributions, (dark grey). The bistable 
parameters are C = 100, 6 = 12, 6c = —0.2, I2 = 40 as in Table II. 

Fig 5. Minimal variance AS'min as a function of the input field squeezing Rtn for p = 
1/2000 [A system] and p = 1/2 [classical 2-level system]. 

Fig 6. Transfer efficiency 1] versus effective cooperativity C for the same systems as in 
Fig. 5. 

Fig 7. Minimal variance versus effective detuning 6 in the case when the ratio rp/70 = 1, 
calculated using the 2-level model with prescription (|6^ and the full 3-level model. 
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Fig 8. Minimal variance versus effective detuning 5 in the case of a closed system [Tp S> 
7o], calculated using the corrected 2- level model and the 3- level model. The small field term 
in Eq. (|63b|) allows for going below the 1 / a/2 minimal variance limit and improving the 
squeezing (up to 35% in this case). 



[1] M. D. Reid, Phys. Rev. A 37, 4792 (1988) 

[2] F. Castelli, L.A. Lugiato, M. Vadacchino, Nuovo Cimento SOC Ital. Fis., B 10, 183 (1998) 

[3] L. Hilico, C. Fabre, S.Reynaud, E. Giacobino, Phys. Rev. A 46, 4397 (1992) 

[4] A. Lambrecht, T. Coudreau, A.M. Steinberg, E. Giacobino, Europhys. Lett. 36, 93 (1996) 

[5] A. Sinatra, J.F. Roch, K. Vigneron, Ph. Grelu, J. Ph. Poizat, K. Wang, P. Grangier, Phys. Rev. 
A 57, 2980 (1998) 

[6] T. Coudreau, L. Vernac, A.Z. Khoury, G. Breitenbach, E. Giacobino, Europhys. Lett. 46, 722 
(1999) 

[7] D.F. Walls and P. Zoller, Phys. Rev. Lett. 47, 709 (1981) 

[8] C.W. Gardiner, Handbook of Stochastic Methods, Springer Verlag (1985); H.J. Carmichael, Phys. 
Rev. A 33, 3262 (1986) 

[9] M. Kitagawa and M. Ueda, Phys. Rev. A 47, 5138 (1993) 

[10] W.M. Itano, J.C. Bergquist, J.J. Bollinger, J.M. Gihigan, D.J. Heinzen, F.L. More, M.G. Raizen, 
D.J. Wineland, Phys. Rev. A 47, 3554 (1993) 

[11] D.J. Wineland, J.J. Bohinger, W.M. Itano, D.J. Heinzen, Phys. Rev. A 50, 67 (1994) 

[12] L. Vernac, M. Pinard, E. Giacobino, Phys. Rev. A, 62, 063812 (2000) 

[13] L. Vernac, M. Pinard, V. Josse, E. Giacobino, Eur. Phys. J. D., 18, 129 (2002) 



25 



[14] A. Kuzmich, K. Moelmer and E.S. Polzik, Phys. Rev. Lett. 79, 4782 (1997) 

[15] J. Hald, J.L. Soerensen, C. Schori, and E.S. Polzik, Phys. Rev. Lett. 83, 1319 (1999) 

[16] L. Vernac, M. Pinard, E. Giacobino, Eur. Phys. J. D, 17,125 (2001) 

[17] C. Cohen- Tannoudji, J. Dupont-Roc, G. Grynberg, Atom-Photons Interactions (Wiley, New- 
York, 1991), p 385 

[18] G.S. Agarwal, R.R. Puri and R.P. Singh, Phys. Rev. A 56, 2249 (1997) 
[19] K. Wodkiewicz and J.H. Eberly, J. Opt. Soc. Am, B2, 458 (1985) 



26 



This figure "figl.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0209023vl 



This figure "fig2.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0209023vl 



This figure "fig3.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0209023vl 



This figure "fig4.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0209023vl 



This figure "fig5.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0209023vl 



This figure "fig6.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0209023vl 



This figure "fig7.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0209023vl 



This figure "fig8.jpg" is available in "jpg" format from: 



http://arXiv.org/ps/quant-ph/0209023vl 



